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Abstract 

An atmospheric transparency was measured using a LIDAR with a pulsed UV laser 
(355nm) at the observation site of Telescope Array in Utah, USA. The measure- 
ment at night for two years in 2007 ^ 2009 revealed that the extinction coefficient 
by aerosol at the ground level is O.OSS^QgJg^^m^^ and the vertical aerosol optical 
depth at 5km above the ground is 0.035!!^og}3. A model of the altitudinal aerosol 
distribution was built based on these measurements for the analysis of atmospheric 
attenuation of the fluorescence light generated by ultra high energy cosmic rays. 

Key words: Ultra High Energy Cosmic Rays, Fluorescence, Atmospheric 
transparency, LIDAR 



1 Introduction 

Ultra High Energy Cosmic Rays (UHECRs) interact with the cosmic mi- 
crowave background and produce a cutoff structure in their energy spectrum 
(GZK cutoff) at the energy around 10^^'^ eV as predicted by Greisen, Zat- 
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sepin and Kuzmin [T][2]. The HiRes group reported that there is a GZK cut- 
off in their observed energy spectrum at the energy it was predicted [3] [4]. 
A similar suppression in the energy spectrum is reported by the Auger ex- 
periment in the Southern Hemisphere [5]. However, the AGASA experiment 
observed a spectrum which continued unabatedly with 11 events beyond the 
GZK cutoff |S]IZ]- To understand this discrepancy, the Telescope Array (TA) 
experiment was constructed in the desert area near the town of Delta in Utah, 
the USA. It is a hybrid detector consisting of Surface Detector (SD) array 
and three Fluorescence Detector (FD) stations. It measures the energy spec- 
trum, anisotropy and composition of UHECRs to identify their origin. The 
TA has three FD stations called "Black Rock Mesa" (BR), "Long Ridge" (LR) 
and "Middle Drum"(MD), which have been installed surrounding the SD ar- 
ray. The BR site (1404 m a.s.l.) is located at the southeast corner of the SD 
array while the LR site (1554 m a.s.l.) is to the southwest. The BR and LR 




Fig. 1. Location of the LIDAR system. A left figure is a map of the TA experiment. 
Little black points are SDs, and three FD stations are shown by BR, LR, and MD. 
A right picture shows the positions of BR-station and LIDAR dome. 
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stations each have 12 telescopes. The MD site (1610 m a.s.l.) is located at the 
northwest, and has 14 telescopes. The FD stations overlook an array of 507 
SDs. 

The UV fluorescence light generated by an air shower is scattered and lost 
along the path of transmission to the telescope. The main scattering pro- 
cesses are Rayleigh scattering by molecules and scattering by aerosols in the 
atmosphere. The Rayleigh scattering process is well understood and the at- 
tenuation length can be calculated using the Rayleigh scattering cross-section 
and the molecular densities of the atmosphere [8] [9] [10]. In order to calculate 
the molecular densities of the atmosphere, radiosonde data from Elko(Nevada) 
are used to obtain temperature, pressure and humidity around the TA obser- 
vatory as a function of height. Sizes, shapes and spatial distribution of aerosols 
around the site are not known, and are variable with time. Therefore, on-site 
monitoring of aerosols is essential in a fluorescence experiment. 

In the TA, we employ a variety of measurements for atmospheric monitoring, 
using two laser systems and a cloud camera. The flrst laser system is the LI- 
DAR (Light Detection And Ranging, we call LIDAR) installed near the BR 
station, which injects a pulsed laser light in the atmosphere and observes the 
back-scattered light at the same location. LIDAR is widely used in ground 
based aerosol measurement. The LIDAR system is operated before the be- 
ginning and after the end of an FD observation, twice a night. The second 
laser system is located at the geographical center of the three FD stations. It 
flres a vertical laser beam of 355nm and the scattered light by the atmosphere 
is observed by each fluorescence detector station. This system is called CLF 
(Central Laser Facility) [II][l2]. We are shooting 300 shots of laser at lOHz 
from the CLF every 30 minutes during FD observation. In addition, we in- 
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(a) LIDAR telescope and laser (b) Block diagram of the LIDAR system 
Fig. 2. The LIDAR system. A left picture is LIDAR's optical system (telescope and 
laser, etc.). A right picture is a connection block diagram of device of LIDAR. 
stalled an infrared CCD camera for cloud monitoring near the LIDAR system, 

and take pictures of the night sky every hour during FD observation [13] . Fur- 
thermore, the weather monitors are installed in each of the FD stations and 
at the CLF to obtain the atmospheric status at the ground level. 

In this paper, we report on the LIDAR system (section [2]), the data set and the 
analysis method (section [3]), the determination of the Rayleigh scattering by 
the radiosonde data (section HI), the result of aerosol scattering by the LIDAR 
data (section [5]), and a model of atmospheric transparency at the TA site 
(section [6]). 

2 LIDAR system 

A photograph of the LIDAR system is shown Fig|2]^a). The LIDAR is com- 
posed of five basic system blocks. Those are (1) a laser, (2) a steerable tele- 
scope, (3) a photo-multiplier tube, (4) a digital oscilloscope, and (5) a control 
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system as shown in Fig|2]^b). We use an air cooled Nd:YAG laser (Orion model 
by New Wave Research) with a third harmonic oscillation module (355nm). 
The laser fires a pulse with the width of 4 ~ 6ns at IHz. The maximum energy 
of the laser pulse is 4mJ. A MEADE LX200GPS-30 telescope with a diameter 
of 305 mm and a focal length of 3048 mm is used with a photo-multiplier tube 
(HAMAMATSU R3479) mounted at the focus of the telescope. A linear range 
of the PMT was checked by simultaneously firing a series of UV-LED pulses 
at the telescope overlaying with the laser shot, and confirming that the over- 
laid signal has the same charge as the LED-only signal without the laser shot. 
The laser is attached to the telescope mount, therefore the direction control of 
the laser-telescope-PMT system can be simply accomplished by commands to 
the telescope mount. The back-scattered photons from the laser are collected 
by the telescope, detected by the PMT, and the signals are digitized with a 
digital oscilloscope (Tektronix 3034B) [11] [13] . In addition, a portion of each 
laser shot is picked off for measurement by an energy monitor. 

By measuring the time structure of the back-scattered photons, we can de- 
termine the atmospheric conditions including aerosol distributions along the 
path of the laser beam. The operation of the LIDAR is composed of horizon- 
tal shots in the north and vertical shots. The vertical operation was made in 
high (~ 4mJ) and low (~ ImJ) energies in order to measure the extinction 
coefficient a over a large range. The horizontal operation uses only the high 
energy. In each vertical and horizontal operation, 500 laser pulses are shot and 
recorded. A total of 1500 shots composed one LIDAR operation. The pedestal 
level is measured in between each shot and is subtracted from the laser shot 
data in order to account for the background. In the following analysis, we use 
an average of the PMT signal profiles in order to reduce shot-to-shot fluctua- 
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Fig. 3. Example of W{t) for horizontal Fig. 4. Example of F{x) as a function of 

shots the distance x. 

tions. 



3 Data set and analysis method 

A waveform W{t) is obtained by averaging the PMT signal for 500 laser shots. 
A range corrected power return F{x) of the LIDAR is defined by the waveform 
W{t) as a function of the distance x from the LIDAR. 

F{x) = W{t)x', x = t^, (1) 

where c is the speed of light and t is the time from the laser shot. The solid 
angle correction is taken into account in the definition of F{x). An example 
of W{t) and F{x) is shown in Figl3]and FigJH 

The extinction coefficient a{x) at a distance x is related with F{x) by the 
LIDAR equation; 



1 dFW _ 1 d/3(.) ^^i^j 



F{x) dx f3{x) dx 

where /3 is the back-scattering coefficient. The factor 2 of a indicates a round 
trip photon propagation. 
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If the atmosphere is spatially uniform along x, such as the case for the hor- 
izontal shots, dp{x)/dx = and F{x) becomes exponential. We can exactly 
solve the equation and obtain the extinction coefficient a{x) (the slope method 

mm)- 

In the case of vertical laser shots, where the atmospheric condition {a and /3) 
changes with x, we need an additional assumption on the relation between a 
and /3 and also a boundary condition to solve the equation. In Klett's method 
[T7][TB], a boundary condition is assumed that a at enough high altitude is 
given by that of Rayleigh scattering because the aerosol scattering can be 
neglected at high altitude. The back-scattering coefficient /3 is assumed to be 
P oc a'^. The value of k is 1.00 for the molecular atmosphere but is known to 
vary significantly (0.67 < n < 1.30) for the atmosphere heavily loaded by the 
aerosol or with rain and snow [H] [19] [20] [21] . However, for a typical desert 
atmosphere of the TA site, we found that the aerosol contribution increases 
towards the lower altitude and that it modifies the value of k from 1.00 at 
the high altitude up to 1.14 near the ground. In solving Eq. (|2]) by Klett's 
method, we assume that (3 is in proportion to a for the whole interval of x. 
We estimate a systematic error caused by this assumption by the iterative 
method (see section |5]). 

The extinction coefficient at the distance x is thus calculated from the acquired 
waveform W{t) by adopting the slope method for horizontal shots or the 
Klett's method for the vertical shots. The obtained extinction coefficient, aobs? 
contains effects of both attenuation processes caused by Rayleigh and aerosol 
scattering. The extinction coefficient by the aerosol scattering a^s is obtained 
using a relation aobs(a;) = aRay(a^) + c^Asi^) where a^ay is the extinction 
coefficient by Rayleigh scattering. 



8 



The LIDAR operation has been performed 515 times from September 2007 
to October 2009. Out of 515 operations, we passed 325 operations to the 
detailed analysis after removing the operations performed under apparent bad 
weather or bad atmospheric conditions, with irregular forms of F{x) caused 
by the existence of clouds and other scatters, and of hardware troubles. Two 
selection criteria on the data quality, the first for the horizontal shots and 
the second for the vertical shots, are then applied on 325 samples. Firstly, we 
discarded the data with Oobs on the ground level larger than O.lkm"^ because 
the atmosphere is too hazy for meaningful measurement. We are left with 
250 operations after this cut. Secondly, we accepted only measurements for 
which the lowest viable data is below 1.0km and the highest viable data is 
above 6.0km. A total of 137 operations have passed the second cut. These two 
conditions indicate a transparent atmosphere and the altitudinal distribution 
of Oobs can be derived from the observed LIDAR data with high accuracy. 
In section [5l we use these 137 LIDAR measurements for characterizing the 
aerosol scattering. 

4 Scattering by atmospheric molecules 

The molecules (nitrogen and oxygen) are the main sources of photon scattering 
in the atmosphere. The scattering processes of photons of wavelengths in the 
range between 300 and 400 nm [22], which is the sensitivity band used in air 
fluorescence experiments, are well described by atmospheric molecules based 
on the Rayleigh scattering theory. We can calculate the attenuation of the 
photon intensity, or the extinction coefficient, from the Rayleigh scattering 
cross-section and the molecular density of the atmosphere. The former can 
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be calculated with an accuracy of 1% [9] and the latter is obtained by the 
measured distribution of atmospheric temperature and pressure employing the 
equation of state [23]. 

We investigated radiosonde data of several stations around the TA site to 
determine the atmospheric conditions. Worldwide balloon-borne radiosonde 
measurements are carried out at more than 900 launch stations twice a day at 
0:00 and 12:00 UTC, and various atmospheric parameters including tempera- 
ture and pressure are recorded along the altitude up to 30 km above the sea 
level. We examined the data from 5 radiosonde stations within 500 km range 
of the TA site, and found that the atmospheric parameters obtained from the 
Elko(Nevada) are most suitable for TA data analysis. Elko is located about 
390km northwest of the TA site. It is at a similar altitude to the TA site and 
has similar climate, temperature, and precipitation. 

In calculation of the extinction coefficients of Rayleigh scattering, we use 
monthly averages of the atmospheric parameters at 00:00 and 12:00. The dis- 
tribution of monthly averages of ftRay from January 2007 to December 2009 
is shown in FigJS] The average of ctRay at the ground level and its standard 
deviation are 0.060 ± 0.002km~^, which translates in an average attenuation 
length of Rayleigh scattering of 16.7km or 1747g/cm^. The systematic error 
of ORay is approximately 1% coming from the inaccuracy of the Rayleigh scat- 
tering cross-section. At the TA site (1400m a.s.L), an average temperature of 
10.2 °C and atmospheric pressure of 863.3hPa at the ground level is measured. 
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Fig. 5. Distribution of monthly averages of aRay(Okm) at A = 355nm. 
5 Scattering by the aerosol component 

The aerosol component consists of large molecules and particles of sizes up 
to ~ 10""^ cm. It influences the light scattering mostly at the low altitude. 
The light scattering by a sphere of a given size is described by Mie scattering 
theory, however, it is difficult to apply the theory for light scattering in the 
lower atmosphere because we do not know shape and the spatial distribution 
of the aerosols a priori. 

The method used in this analysis is as follows. First, a total extinction co- 
efficient Oobs that is the sum of Rayleigh scattering and aerosol scattering is 
determined using the LIDAR data, in a unit of km~^, as a function of the 
height. Next, we calculate the extinction coefficient of Rayleigh scattering 
ctRay using the radiosonde data. Finally, the extinction coefficient of aerosol 
scattering a as is calculated by assuming aohs{x) = aRa.y{x) + aAs{x). 

5.1 Atmospheric transparency in the vicinity of the ground surface 

The aerosols are more prevalent at the ground surface and decrease as a func- 
tion of the height. They are influenced easily by the wind, ground dryness and 



11 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 O.OS 0.09 



Fig. 7. Distribution of the maximum 
Fig. 6. Distribution of the Q;As(Okm). height below which the aerosol contribu- 
tion is measurable. 

other local factors. Figure [6] shows the distribution of oas at the ground level 
obtained from horizontal measurements by the slope method. The median of 
OAS (0km) and its range (Icr) of distribution over 137 LIDAR operations are 
0.033lQ;oi2kiii~^. The attenuation length of aerosol scattering at the ground 
level was determined to be 30lxokni. The systematic error of oas by the lin- 
earity calibration of the PMT gain is 1*17%. 

5.2 Atmospheric transparency of each altitude 

The altitudinal distribution of Oobs per 250m is calculated by Klett's method 
from the data of a vertical shooting. Examples of extinction coefficient ob- 
tained from the measurement for the different distributions of aerosol scat- 
tering are shown in Fig|8](a)~(d). The Oobs obtained at the ground level by 
the horizontal shooting and the distribution of Oobs measured by the vertical 
shooting are shown in the closed circle and the crosses in Fig|8] respectively. 
For vertical shots (crosses), we used data only above 750m for which the lin- 
ear response of the PMT was confirmed by the UV-LED calibration. The data 
points above 8km are removed because the collected photon signal is too small. 
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(c) Aerosol distributed up to high altiutde (d) Aerosol distributed at both altitudes 
Fig. 8. Extinction coefficients as a function of the height from the ground level. 
The solid line of Fig|8] represents aRay calculated for pure Rayleigh scattering 
(see section Hj). Finally the open square boxes in Fig J8] represent the altitudinal 
distribution of aAs(a;)- 



The systematic error of aAs{x) is from the linearity calibration of the 
PMT, lg% from the assumption of k = 1 for Klett's method. The latter is 
estimated by iteration; at first, the aerosol contribution was obtained by the 
Klett's integration with the assumption /3 oc a'^^^, then the relation between 
a and /3 was recalculated at each altitude assuming the aerosol phase function 
obtained by the HiRes experiment [23] and the second Klett's integration was 
performed. We found the second integration gives a correction of +5% or less 
at any altitude of the typical atmosphere. 
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Fig. 9. Median of TAs{h)- Fig. 10. Distribution of tas at 5.0km. 

In order to understand the aerosol contribution with respect to the height, 
we define a height below which the a as is greater than O.Olkm"^, i.e. aerosol 
scattering is no longer negligible. The distribution of such height is shown in 
FiglTl Although the distribution reaches as high as 5km, it is less than 2 km 
in most cases (71%). We found, for 2% of the cases, the contribution from the 
aerosol scattering was completely negligible. 



5.3 Vertical Aerosols Optical Depth 



In this section, we evaluate a transmittance of the atmosphere in terms of 
the Vertical Aerosols Optical Depth (VAOD). The VAOD at the height of h, 
TAs{h), is defined by the integration of a^s from to /i as follows, 

TAs{h)= raAs{h')dh'. (3) 
Jo 

The attenuation factor for photons which vertically propagate in the atmo- 
sphere until height h is given by exp {—TAs{h)} . The TAs{h) was calculated as 
a function of height using aAs{x) night by night. 

In the integration of Eq. ([3]), the values of aAs(^) between h = and hmin are 
assumed to be the average of two curves; 
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1. Exponential curve connecting aAs(O) and oas (^mm) 

2. Extrapolation of the exponential curve fitted with the aAs(^) between hmin 
and hmin + 1km 

where hmin is defined as the minimum height for which oas can be determined 
by vertical shots. These two curves are shown by two broken lines in FiglHl 

The median of TAs{h) is shown in FigJHl The distribution of tas at 5km above 
the ground is shown in FigJTOl The median of tas at 5km and its range (Icr) 
of the distribution are 0.035lo;oi3. The systematic error of the VAOD was 
evaluated by taking these two curves as the upper and lower limits, and is 
estimated to be ±6%. Additional systematic errors are 11%% from the linearity 
correction of the PMT response and to% by the assumption of k = 1.0 in the 
Klett's integration method. The total systematic error of VAOD is ^17% added 
in quadrature. 

The seasonal variation of tas distribution is shown in FigJTTl The data of 
winter and summer include the data for five months that centered January 
and July, respectively. The mean values of tas and their distribution range 
(1 a) are O.OSgto;"?" in summer and 0.025lo:oo? in winter. The effect of the 
aerosol component in summer is 1.5 times greater than that in winter. 

6 Model of transparency distribution 

In this section, a model of the altitudinal distribution of the atmospheric 
transparency based on the aAs(^) measurement is presented. Figure [12] shows 
the distribution of oas averaged over 137 LIDAR operations as a function of 
the height h above the ground level. The median of ctAs shown in Fig{T2] can 
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(a) Median of TAs{h) for different seasons (b) Distribution of tas at 5km above 

ground level for summer (above) and winter 
(below) 

Fig. 11. Aerosol scattering status for different seasons, 
be well expressed by two- component exponential function; 



aAs{h) = 0.017exp(-V2.1) +0.016 exp(-/i/0.1) 



(4) 



where h is the height from the ground level of LIDAR observation (1400m 
a.s.l.) in the unit of km. The first term of the function (jl]) is responsible for 
the distribution of the aerosol at high altitude {h > 0.75km) whereas the lower 
altitude distribution is expressed by the super position of two terms. Especially 
at the ground level (/i ~ 0), two terms have nearly equal contributions. The 
functional from of aAs{h) is plotted in Figll2]as the solid line. The integral of 
C(As{h) is overlaid as a solid line in Figj9]and it reproduces the median value 
of measured tasW w^U. 



7 Conclusion 

The LIDAR atmospheric monitoring system of TA began operation in Septem- 
ber 2007 and continues until today. Analysing the data from the LIDAR, we 
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(a) linear scale (b) log scale 

Fig. 12. Distribution of oas as a function of the height from the ground level. At 
each height, crosses show a median of oaS; broken error bars indicate the range 
(68%) of its distribution. 



have found that the extinction coefficient by aerosol at the ground level is 
OAS (0km) = 0.033lo;oi2kni~^ corresponding to SOligkm in the attenuation 
length where the error (±) corresponds to the range of distribution (Icr). The 
maximum altitude at which the contribution of the aerosol scattering was ob- 
served, excluding the cloud, is 5km. We obtained the VAOD at each altitude 
/t[km], integrating aAs(^) from to h. The median of the VAOD at 5 km 
above the ground is rAs(5km) = 0.035lo;o|3 with a systematic error of li7%. 
The effect of the aerosol in summer is found to be 1.5 times greater than that 
in winter at the observation site. The average altitudinal distribution of oas 
above the ground level (1400m a.s.l.) is well reproduce by the super position 
of two exponential functions; 



aAs{h) = 0.017 exp(-/i/2.1) + 0.016 exp(-/i/0.1). 
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